Nonequilibrium phase transition in a model for social influence. 
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We present extensive numerical simulations of the Axelrod's model for social influence, aimed at 
understanding the formation of cultural domains. This is a nonequilibrium model with short range 
interactions and a remarkably rich dynamical behavior. We study the phase diagram of the model 
and uncover a nonequilibrium phase transition separating an ordered (culturally polarized) phase 
from a disordered (culturally fragmented) one. The nature of the phase transition can be continuous 
or discontinuous depending on the model parameters. At the transition, the size of cultural regions 
is power-law distributed. 
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Recently, the study of complex systems entered social 
science in order to understand how self-organization, co- 
operative effects and adaptation arise in social systems 
jij. In this context the use of simple automata or dy- 
namical models often elucidates the mechanisms at the 
basis of the observed complex behaviors (l]J^] . 

In this spirit, R. Axelrod has recently proposed an in- 
teresting model to mimic how dissemination of culture 
works HI]]. Culture is used here to indicate the set of 
individual attributes, such as language, art or technical 
standards subject to social influence, i.e. that can be 
changed as effect of mutual interactions. The automa- 
ton does not consider the effect of central institutions or 
mass media and focuses on the self-organization result- 
ing from a simple local dynamics representing the social 
influence. This dynamics is assumed to satisfy two sim- 
ple properties: i) individuals are more likely to interact 
with others who already share many of their cultural at- 
tributes; ii) interaction increases the number of features 
that individuals share. Starting from an initial state with 
features distributed randomly this leads to the formation 
and coarsening of regions of shared culture. 

In this Letter we carry out an accurate numerical anal- 
ysis of Axelrod's model that unravels a remarkably rich 
behavior, not detected in previous investigations. De- 
pending on the initial degree of disorder, the model un- 
dergoes a phase transition separating an ordered from a 
disordered phase. The ordered phase is characterized by 
the growth of a dominant cultural region spanning a large 
fraction of the whole system. On the contrary, in the dis- 
ordered phase the system freezes in a highly fragmented 
state with a nontrivial distribution of the sizes of cultural 
regions. Such a fragmented configuration is reached in a 
finite time, which diverges at the phase transition. In the 
whole ordered phase instead, the coarsening process lasts 
for a time proportional to the system size, before freezing 
into the culturally polarized state. Interestingly, the na- 
ture of the transition turns from continuous to discontin- 
uous when the number of cultural features is increased. 
Close to the transition, the distribution of region sizes 



follows a power law in agreement with the results of Ref. 
Pj on the statistical analysis of the diversity of languages. 
Some of these features are captured by a mean field ap- 
proximation that we discuss below. 

Axelrod's model is defined on a square lattice of linear 
size L. On each site i there is a set of F integer variables 
aij which define the cultural "features" of the individu- 
als living on that site. In the original model, each feature 
/ = 1, . . . , F on each site i is initially drawn randomly 
from a uniform distribution on the integers between 1 
and q. The parameter q is a measure of the initial cul- 
tural variability (i.e. disorder) in the system. Here, we 
relax the constraint of integer q by using a Poisson dis- 
tribution Prob((7ij = k) = q k e~ q /kl, that allows q to be 
a positive real number. Though the results are qualita- 
tively the same, our choice is more convenient to study 
the behavior of the model as q varies. 

At each time step, a pair of nearest neighbor sites i 
and j is randomly chosen. A feature / is chosen ran- 
domly and if cr^j ^ (Tjj nothing happens. If instead 
&i.f = Cjj then an additional feature /' is randomly 
chosen among those taking different values across the 
bond, Uijt ^ o~jj'. Such a feature is then set equal: 
Cj,/' —* o~\ f = ■ A "sweep" of the lattice, i.e. I? 
such time steps, defines the time unit as usual. Axelrod's 
model (at least in the original formulation) can be seen 
as F coupled voter models ||. 

During the dynamical evolution the total diversity, 
measured as the number of different values of a fea- 
ture / which are present in the system, always decreases. 
Clearly if all features are equal across a bond {o~ij = Oj,/, 
V/) or if they are all different (<Tij ^ <jjj V/) no change 
can occur on the bond ij. A configuration such that, 
on each bond, either all features are equal or they are 
all different is an absorbing state: Dynamics will stop if 
such a state is reached. There are clearly many absorb- 
ing states: The dynamics on any finite lattice converges 
to one of them. The final state can be characterized by 
the distribution of cultural region sizes, where a region is 
defined as the connected set of sites sharing exactly the 
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same features. 

The dynamical evolution is characterized by the com- 
petition between the disorder of the initial configuration 
and the ordering drive due to the local social interactions. 
It is intuitively clear that when q is small the initial state 
is almost completely uniform, whereas for large values of 
q almost all sites have features atj totally different from 
those of their nearest neighbors. In the two cases we 
expect the system to converge to a uniform or a highly 
fragmented state, in which interaction or disorder dom- 
inate, respectively. In order to understand how these 
two limit situations are connected as q varies, we have 
simulated the Axelrod's model for a number of features 
ranging from F — 2 to 10 and sizes up to L = 150. We 
first discuss the dependence of the final state on the pa- 
rameters q and F and then the dynamical behavior of 
the model. 

The frozen state. In any finite lattice the dynamics 
converges to a frozen absorbing state. The existence of a 
transition in the properties of the final absorbing states 
is very clear from the plot (Fig. [l]) of the average size of 
the largest region (s max ) as a function of q for F = 10: 
For q = q c ~ 300 we observe a sharp transition char- 
acterized by a sudden drop of (s mayi )/L 2 which becomes 
steeper and steeper for increasing sizes L. This points to 
the existence of a transition between a "culturally polar- 
ized phase" for q < q c , where one of the regions has a 
size of the order of the whole system, and a "culturally 
fragmented" phase, where all domains are finite. The 
transition is of the first order, with the size of the largest 
region having a finite discontinuity at q = q c . Note also 
that, for q < q c the largest domain approaches a unitary 
density, i.e. it invades the whole system. This scenario 
holds for all values F > 2 investigated. 

The situation is different for F — 2 (inset of Fig. II]): 
The fraction occupied by the largest cluster (s max )/L 2 
vanishes continuously as q — > q~ . The difference be- 
tween F = 2 and F > 2 is confirmed and clarified by the 
study of the size distribution of cultural regions at the 
transition. Let Pl(s, q) be the probability distribution of 
regions of s sites in a system of size L. The cumulated 
distribution is plotted in Fig. || for several F and values 
of q around the transition. Fig. || shows that Pi,(s,q) 
decays as a power law s~ T and that the exponent r is 
universal (r ~ 2.6) for F > 2 but takes a different value 
for F = 2 (r ~ 1.6). In particular we find Tp>2 > 2 and 

T F =2 < 2. 

The different nature of the transition for F equal or 
greater than 2 can be related to the exponent r. Let 
N(q, L) be the total number of regions in the system. 
Requiring the total area to be L 2 leads to 

oo 

L 2 = N(q, L)(s) = N(q, L) £ a P L (s, q). (1) 

s=l 

For q > q c there are N(q, L) ~ L 2 domains of finite size 
and the sum on s is finite as L — > oo. On the other hand, 
for q < q c , there are few small domains and a large one 



of size s max ~ L 2 . Hence the probability distribution 
can be written in the generic scaling form P^s^q) — 
s~ T J-(s I Sco) + A(q)d s , Sin ^ where s co is a cut-off scale, 
the function !F(x) is constant for x <C 1 and decays very 
rapidly for x S> 1 and A(q) — for q > q c . 

The divergence L 2 in the l.h.s. of Eq. (|) is matched, 
for q < q c , by the component j4<5 S)Smax in P^(s,q), with 
•Smax ~ L 2 . The nature of the transition is identified by 
the behavior of A(q) for q — > q~ . For r < 2, similarly 
to what happens in percolation theory , the transition 
occurs through the divergence of s co and hence of a cor- 
relation length, as q — > q+. This causes the divergence 
of (s) in Eq. (Q) (because sPi(s,g) ~ s 1_r with r < 2), 
which matches the divergence of the l.h.s. of Eq. (0) as 
L — > oo, through usual finite size scaling arguments [0. 
Indeed, at q = q c the cut-off diverges with the system size 
as s co ~ L D . On the other side of the transition a similar 
divergence of the cutoff s co occurs and the amplitude A 
must vanish as q — > q~ . This scenario is typical of second 
order phase transitions and agrees with the value r < 2 
found for F — 2 and Fig. @. 

When r > 2 the previous scenario cannot hold: For 
q > q c , even if s co diverges, the sum on s in Eq. (0) 
remains finite. Hence as q — > g+ the number of domains 
N(q, L) must remain of order L 2 . On the other side of the 
transition, the L 2 term in the l.h.s. of Eq. (0) can only 
be matched by the term s max ~ L 2 . Since no divergence 
arises from the sum on small components as q — > q~ , the 
amplitude A of the large component must remain of order 
one in this limit. We then conclude that, as q crosses q c , 
the nature of the distribution changes abruptly for r > 2. 
In particular the amplitude A(q) exhibits a discontinuous 
jump across the transition [g. 

Dynamics: In order to investigate the model dynamics 
we study the density n a of active bonds. An active bond 
is a bond across which at least one feature is different and 
at least one is equal, so that there can be some dynamics. 
This quantity is indicative of the dynamical state of the 
system, being zero in frozen configurations. In Fig. [| we 
show the behavior of function of time for different 

values of q and F — 10. For large values of q, after a 
short initial transient, the density of active bonds decays 
rapidly and the system locks into a frozen configuration 
in a finite time. For small q instead, n a (t) displays a slow 
decay with a large majority of bonds in the active state 
until it falls abruptly for long times. As the system size 
is increased, the slow decay extends to longer and longer 
times: The cut-off time, at which activity suddenly dies 
scales as t co ~ L 2 as shown in the inset of Fig. ||. For 
a large range of intermediate q values, n a first decreases 
almost to zero but then rises again towards a peak of 
activity, from which the slow decay begins. 

The dynamics is essentially a coarsening process of ho- 
mogeneous regions. When the process lasts only for a 
finite time, as for q > q c , it gives rise to regions of finite 
size. On the contrary, for q < q c the coarsening process 
goes on for a time t co ~ L 2 and produces a region of 
size comparable with the whole system. The exponent 
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z = 2 relating t co to L is the same found in noncon- 
served phase-ordering || and in particular in the two- 
dimensional voter model ]To| , which also exhibits a slow 
decay of the density of active bonds jll| . 

For an infinite system (L — > oc) and q < q c the sys- 
tem is indefinitely in a coarsening state. Therefore we 
can also define the transition as separating two different 
dynamical regimes, with a fast decay of n a lasting for a 
finite characteristic time for q > q c and a slow, infinitely 
long decay for q < q c . For F = 2 the behavior is qualita- 
tively the same: Noticeably, no evident signature of the 
different nature of the transition (continuous vs. discon- 
tinuous) can be inferred from the dynamical evolution. 

The dynamical behavior of the model can be studied 
within a single bond mean field treatment. Let P m {t) be 
the probability that a randomly picked bond is of type m 
at time t, i. e. m features across the bond are equal and 
F — m are different. At t = 0, since features are assigned 
uncorrelated random values, we have P m (0) = ( ) p™{l — 
Po) F ~ m where po = Prob[<Jij = (Tjj] is the probability 
that two sites have feature / with the same value at t = 
In the mean field approximation, P m satisfies the master 
equation 

k=l 



(g (PnW r [% - P m WH 
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where g is the lattice coordination number and W, 
is the transition probability from a n-type bond to a m- 
type bond due to the updating of a fc-type neighbor bond. 
This equation describes how the number of bonds of type 
to is affected by the dynamics: k/F is the probability 
to select a bond of type k and one of the k features which 
are equal across it. If k = m — I a new bond of type m is 
created whereas if k = m one bond of type m is deleted. 
This explains the first two terms on the r. h. s. of 
Eq. (|J) . But the change of a feature of the site may also 
affect the state of the other g — 1 bonds connecting such 
site to its neighbors. The remaining terms in Eq. (j^) 
take into account this possibility. We can compute the 

(k) 

probabilities Wn,m by analyzing in detail each possible 
process. Let us consider e.g. a type k = 1 bond adjacent 
to a type n = one. Without loss of generality we can 
let the features on the extreme sites of the 1 bond be 
(0, 0, ... , 0) and (0, 02, . . . , (Jf), with <Tj ^ 0. The latter 
site is shared also by the bond whose other site has 
features (a[, (r' 2 , . . . , a' F ) with o' x ^ and a'j ^ <Tj for 
j > 1. When the 1 bond becomes a type 2 bond, o 2 — > 
the neighbor bond can either remain a bond, if a' 2 
or it can become a 1 bond if a' 2 = 0- In the spirit of the 
mean field approximation we introduce the probability 
p = Prob(<72 = 0) which now becomes time dependent, 



as we shall see. Then W^l = 1 — p and = p. In 

much the same way we can compute the other transition 
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rates. For example, if F = 2, the only non zero elements 
are 



W, 
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The system of equations (g) is closed once the dynam- 
ics of p(t) is specified. The simplest such equation, in 
the mean-field spirit, is p = kP^/F. This amounts 
to say that between any two sites there is a bond, and 
hence that the probability that a feature across that bond 
takes the same value can be expressed in terms of the P m 
p2[ . The numerical integration of the mean field equa- 
tions yields the phase diagram of the model as shown 
in Fig. |^, which exhibits a phase transition. The order 
parameter 77 a = X^fcL] 1 Pk undergoes a discontinuity at 
q = q c , as it jumps from a finite value for q < q c to zero 
for q > q c . The mean field fairly reproduces also the dy- 
namical evolution of the model as can be seen from the 
inset of Fig. |[ 

Despite the extreme simplicity of Axelrod's model, it 
is tempting to compare the picture derived in this Let- 
ter with some recent analysis on the diversity of lan- 
guages ||, which to some extent can be considered as 
an indicator of cultural homogeneity. Gomes and coau- 
thors, by analyzing the statistics of more than 6700 lan- 
guages, have found that the number C of linguistically 
homogeneous regions inhabited by N individuals scales 
as C ~ N~ T , with r ~ 1.5 for populations smaller than 
6 x 10 6 individuals. For larger populations the expo- 
nent changes suggesting a different dynamical evolution. 
These findings can be compatible with our continuous 
phase transition scenario (r ~ 1.6) assuming the equiva- 
lence between language and cultural features. However it 
remains unclear why language spreading is self-organized 
close to the transition point or why the value F = 2 
should be relevant for the process. 

In summary, we have presented numerical simulations 
of the Axelrod's model for social influence. Social inter- 
action that tends to make culture homogeneous competes 
with the disorder introduced by the number of different 
traits each cultural feature initially has. We find a critical 
value separating two phases in which one of the above el- 
ements (interaction or disorder) respectively dominates. 
The nature of the phase transition and the emergent col- 
lective behavior is analyzed by common tools of statisti- 
cal physics. We find that the transition is continuous or 
discontinuous depending on the model parameters. Our 
study shows that the use of concepts and methods devel- 
oped in physics may be of help in the context of social sci- 
ences. For instance, the prevalence, for small values of q, 
of one of the initially equivalent cultures can be regarded 
as a spontaneous symmetry breaking due to stochastic 
fluctuations. It would be interesting, in the future, to 
analyze extensions of the simple Axelrod's model to take 
into account additional ingredients like the presence of 
migrations or the effect of geographical barriers. 
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FIG. 1. Behavior of s m ax/£ 2 vs q for three different system 
sizes and F = 10. In the inset the same quantity is reported 
for F = 2. 
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FIG. 2. Cumulated distribution of region sizes for q m q c , 
L — 100 and several values of F. 
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FIG. 3. n a (t) for F = 10, q = 1, 100, 200, 230, 240, 250, 
270, 300, 320, 400, 500, 10000 (top to bottom), L = 150. The 
inset reports the dependence of the freezing time t co on L for 
F = 10 and g = 100 < q c . The bold line has slope 2. 
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FIG. 4. Phase diagram within the mean field approxima- 
tion for F = 3. In the inset the mean field behavior of n a (t) 
for q = 15 < q c (dotted), q = 20 ~ q c (full) and q — 25 > q c 
(dashed). In the mean field description the coarsening phase 
lasts indefinitely because it is implicitly assumed an infinite 
size system. 
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